Evidence for topological nonequilibrium in magnetic configurations 
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We use direct numerical simulations to study the evolution, or relaxation, of magnetic 
configurations to an equilibrium state. We use the full single-fluid equations of motion for 
a magnetized, non-resistive, but viscous fluid; and a Lagrangian approach is used to obtain 
exact solutions for the magnetic field. As a result, the topology of the magnetic field remains 
unchanged, which makes it possible to study the case of topological nonequilibrium. We find 
two cases for which such nonequilibrium appears, indicating that these configurations may 
develop singular current sheets. 
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I. INTRODUCTION. 

Formation of singularities, or current sheets, is one of 
the striking features of astrophysical as well as tokamak 
plasmas [jj]. Such singularities are key to understanding 
active phenomena related to fast magnetic field recon- 
nection 0], f|. For example, fast dynamos rely on fast 
reconnection of magnetic field lines , || . Despite their 
importance, key issues related to current sheet formation 
are still not well understood. Supposing, e.g., that they 
are formed due to instabilities, one has to assume that 
fluid dynamical processes are able to slowly deform equi- 
librium magnetic field configurations (and thereby build 
up regions of field gradients) without significant recon- 
nection until a marginal state is reached. At this thresh- 
old, instability-driven reconnection would then lead to 
release of the stored free energy on the (observed) time 
scales thought to be too short to be consistent with, for 
example, Sweet-Parker reconnection ||, m. However, it 
has been long recognized || that in the presence of recon- 
nection, it is not obvious how one can attain (meta)stable 
configurations which store significant free energy. Fur- 
thermore, it is not clear why reconnection would not 
simply return the system to the marginal state, thus re- 
leasing only a small fraction of the available free energy. 

In this paper, we explore one possible solution to these 
puzzles: We consider specific magnetic field configura- 
tions which could arise from a slow evolution of (sta- 
ble) quasi-equilibria, and then examine their subsequent 
(unforced) evolution. Our aim is to show that there ex- 
ist configurations that evolve initially on the slow rate, 
but that can reach a point at which spontaneous current 
sheet formation occurs. These configurations have been 
referred to as "topological nonequilibria" (TN) M, 0, 
and lead to situations in which the topology of the field 
is such that in a relaxed equilibrium state it inevitably 
contains discontinuities. TN results in spontaneous re- 
connection, because no external forces are involved; and 
in the cases we shall examine, the result is that extraction 



of all of the available free energy becomes possible. 

Finally, we note that an important aspect of this prob- 
lem relates to the fact that there is a direct correspon- 
dence between magnetostatic equilibria and steady Eu- 
ler flows, as pointed out by Moffatt ||]; this problem is 
therefore closely connected to the possible formation of 
singularities in hydrodynamics; see also M, p0[. 



II. DESCRIPTION OF THE APPROACH. 
A. The idea of topological nonequilibrium. 

The main ideas of topological nonequilibrium (hence- 
forth, TN) were formulated rigorously by Moffatt ||, M, 
]To| | . Consider an ideally conductive viscous (ICV) flow. 
We restrict ourselves to incompressible flows. Starting 
with initial magnetic field B(x, t = 0) = Bo(x) of ar- 
bitrary topology, one expects that such a configuration 
will relax to a static state, with zero velocity field, and 
nontrivial magnetic field B^. The latter configuration is 
then called 'topologically accessible' because the field's 
topology does not change during this frozen-in evolu- 
tion. If this relaxed equilibrium state contains disconti- 
nuities, then all of the states in the evolution are referred 
to as TN. It may be expected in a realistic situation, 
when small but finite resistivity rj is taken into account, 
that these discontinuities evolve into finite-width current 
sheets, resulting in efficient reconnection and dissipation 
of the magnetic field. Unfortunately, there are only a 
few special cases for which it is possible to demonstrate 
that TN exists In this paper, we restrict ourselves 
to analysis of two-dimensional configurations, and study 
the evolution of two generic field configurations which 
can lead to TN. 

Of course, in general, there is no reason that a given 
initial configuration is at equilibrium. However, one 
would normally expect that, after relaxation, such a con- 
figuration would evolve to attain a smooth equilibrium, 
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and that the magnetic field evolution subsequently stops. 
Some initial field topologies, however, cannot possibly re- 
lax to a smooth equilibrium, resulting in TN. It is obvious 
that use of the word "nonequilibrium" is not strictly cor- 
rect, because in the final state the field is at equilibrium 
as long as the diffusivity vanishes exactly. However, in 
the spirit of maintaining already existing tradition, wc 
retain this terminology. 

As a result of relaxation, the magnetic field will reach 
an equilibrium state. In two dimensions, B x = d y A, 



By 

equation 



-d x A, and the flux A obeys in equilibrium the 
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dP(A) 
dA ' 



(1) 



where P = 
see, e.g., p 



p + B 2 J(8n), p = p(A), and B z = B Z (A); 
]. As an aside, we note that if the pressure p 
can be neglected, then this equilibrium is force-free; this 
may occur in specific applications such as in the solar 
corona. In addition, the total pressure P (Eq. (|l|)) ought 
to substantially exceed the transverse magnetic energy 
B 2 + By in order to justify the incompressibility assump- 
tion for the evolution to this equilibrium state. 

Equation (Q) is trivially satisfied in the one- 
dimensional case. To start with, suppose that A is a 
function of x only, corresponding to straight field lines 
parallel to y axis. An initial arbitrary distribution is 
generally not at equilibrium. However, after relaxation, 
the field reaches the well-known equilibrium 
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+ P(x) = const, 



(2) 



automatically satisfying ([!]). The same is true for ax- 
isymmetric configurations, when A — A(r), and the field 
consists of concentric circles. 



FIG. 1. Sketch of the type A topology configuration: closed 
nested magnetic field lines. The dashed line here, as in all 
the figures below corresponds to a field line with vanishing 
magnetic field strength. 

Going to two dimensions complicates the problem con- 
siderably. Consider first a configuration with closed 
nested field lines, so that A(x, y) has one maximum (min- 
imum). The configuration is depicted in Fig. 1, and we 
will refer to it below as "case A". Of course, the ax- 
isymmetric configuration is topologically accessible from 
this configuration, and therefore it can reach equilibrium. 
The question, however, is if this equilibrium is unique. 



-0.5 0.0 0.5 





-1.5 -1.0 -0.5 0.0 0.5 



FIG. 2. Ellipse-shaped configuration placed between two 
horizontal magnetic walls; the initial configuration was gen- 
erated by considering two families of (parametric) curves, i.e., 
straight lines and ellipses, and joining them as shown in panel 
(a). As a result of evolution of the configuration shown in 
panel (a), the field evolves such that it is pushed to the walls 
to form discontinuities, as shown in panel (b). Dashed lines 
correspond to Bx = 0. Both panels present results of numer- 
ical simulations, but only selected field lines are depicted for 
illustrative purposes. 

If there exists a magnetostatic equilibrium with type 
A field topology with essentially arbitrary field line ge- 
ometry, e.g., with elliptic field lines, as in Fig. 1, then 
we would expect an arbitrary type A configuration to 
relax to this equilibrium without dramatic changes in 
its geometry. However, suppose this equilibrium exists 
only in axisymmetric form, i.e., can only be realized with 
concentric (field line) circles; then, if this configuration 
were placed between magnetic "walls" (such as regions of 
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strong magnetic fields in the solar corona, or solar wind) , 
as in Fig. 2a, we would expect the formation of discon- 
tinuities (because it would not be possible to evolve to 
the equilibrium state). Furthermore, if we allowed for 
nonvanishing diffusion, then such a configuration would 
not settle down until all magnetic lines are reconnected, 
and the bubble seen in Fig. 2 disappears entirely. 

It is useful to expand slightly on the astrophysical rel- 
evance of this case. Our point is that "case A" shown in 
Fig. 2a can be regarded as an abstraction of a com- 
monly expected field configuration in the solar atmo- 
sphere: Consider the emergence of a magnetic flux tube 
from the solar interior to the corona, where it enters 
a highly conducting medium already suffused by pre- 
existing magnetic fields. If one abstracts such an emerg- 
ing flux tube as a rising cylinder, then the expected field 
topology in planes perpendicular to the tube axis should 
be similar to case A: The nested closed field lines in such 
planes then represent the toroidal field component of the 
emerging flux tube; and the magnetic "walls" shown in 
Fig. 2a represent the projections in such planes of the 
magnetic fields of the surrounding magnetized coronal 
plasma. 
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FIG. 3. Simulations for the (initially continuous) rosette 
structure, as depicted in panel (a). As a result of the re- 
laxation, this configuration evolves into a field containing a 
discontinuity between the two magnetic islands, as shown in 
panel (b) . The presence of the external zero line (dashed line) 
is vital for TN. 

The second type of field topology we consider below 
is what we call type B; this more complicated topol- 



ogy is a "rosette structure" (Fig. 3a), which has been 
investigated experimentally p3[ . In terms of the flux 
function A(x,y), this configuration consists of two max- 
ima, e.g., two "mountains", surrounded by a pedestal, 
i.e., two magnetic islands surrounded by closed magnetic 
field lines going around the two islands; the field van- 
ishes outside the zero-line. If the type B topology cannot 
exist in smooth equilibrium, then a current sheet devel- 
ops, resulting in efficient rcconncction of field lines until 
all field lines of the islands are reconnected, and even- 
tually only one island remains, of the topology of the 
type A. In contrast, if this kind of topology does ex- 
ist in smooth equilibrium, then nothing dramatic would 
happen, and the configuration would relax to this equi- 
librium without any discontinuities. The astrophysical 
context in which this type of configuration may be cre- 
ated is similar to that just described above: consider the 
emergence of two adjacent twisted solar flux tubes into 
a non-magnetized ambient corona; again, the field struc- 
ture in a cross-section perpendicular to the tube axes will 
appear as shown in Fig. 3a. Thus, in both cases A and 
B, we are dealing with the generic case of bounded mag- 
netic flux systems (i.e., systems of magnetic field lines 
which lie within a finite bounding surface on which the 
field vanishes) , which can be regarded as abstractions of, 
for example, isolated flux tubes emerging into the highly 
conductive solar corona. 

Generally, finding TN states is far from trivial. To 
illustrate, let us return to the type A topology. The ax- 
isymmetric equilibrium solution is not unique. For ex- 
ample, one can construct a solution to ([!]), 

A(x, y) = sin kx sin ky, 

depicted in Fig. 4a, which has the same topology of field 
lines as depicted in Fig. 1. This asymmetric field is at 
equilibrium, so that the general answer to the question of 
whether, say, elliptic configurations of the form shown in 
Fig. 1 can be at equilibrium, is affirmative. An analogous 
construction can be carried our for the topology of type 
B; in this case, the solution of (nl) can be constructed as 



A, 



^A n e ik ^ )x+ik ¥ )v , 



where {k x n) f + (k y n) ) 2 = constant (see, e.g., @, §). 
This example is depicted in Fig. 4b; the rosette structure 
shown is at equilibrium without any discontinuities. 

The situation changes if a zero-line (a line where = 
0, Bj_ = {B x , By}, and generally B z ^ 0) is present, such 
as the dashed line shown in Fig. 1 for the type A, and 
in Fig. 3 for the type B field topology. The zero-line 
possess two remarkable properties. First, the magnetic 
field remains zero on this line in the presence of ICV 
flow. Thus, if we write the ideal induction equation in 
the form, 



dB 

~dt 



(B • V)v 
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which in 2D reads, 
dB 



dt 



(B_l • V)v, 



dB* 
dt 



= 0, 



(3) 



it is easy to see that because the left-hand side describes 
transport of any fluid element (in particular, of the zero- 
line) by the motion, and because the right-hand side cor- 
responds to change of the field along the Lagrangian tra- 
jectory (and as the right-hand side vanishes on the zero- 
line), this equation will preserve the property = on 
the zero-line. 




- (°) 



0.0 0.2 0.4 0.6 




FIG. 4. Two examples of field configurations which share 
the field topology of the two cases (A and B) we are studying, 
which do not have any simple symmetries, but are neverthe- 
less smooth equilibria, (a) Type A (marked by the thick line); 
(b) Type B (rosette structure, marked by a thick line). 

Second, if the zero-line has a constant (along the line) 
curvature, e.g., is a straight line or a circle, and if the 
field is also analytical, then the entire configuration will 
have the same geometry as the zero-line. In other words, 
if the zero-line is a straight line, then all other field lines 
are straight as well; alternatively, if the zero-line is a 
circle, then the analytical equilibrium configuration con- 
sists of concentric circles. The proof is easily constructed 
by expanding A(x, y) in the vicinity of the zero-line [ fT2"| . 
Note, however, that the constant curvature zero-line is a 
special case (although it can be regarded as a represen- 
tation of the emergence of magnetic flux on, for example, 
the solar surface, in which geometrically symmetric flux 
bundles straddle the separating "neutral line"); in gen- 
eral, the zero- line is arbitrary in shape, as shown in Figs. 
1-3. Nevertheless, we may conjecture that the zero- 
line imposes a severe constraint on the geometry: That 
is, we conjecture that the existence of this line results in 
unique (smooth) solutions of the equation ([!]) in the form 
of magnetic field lines with constant curvature [fl2|| . 

One of the considerations in favor of this conjecture 
is as follows. Without loss of generality, A = out- 
side the configuration, and thus A = on the boundary 
(whose shape is as yet unspecified) , corresponding to the 
Dirichlet problem for equation (fy). On the other hand, 
because B x = B y = on the same boundary, we have 
d n A — 0, corresponding to a Neumann problem. The 
problem is thus over-constrained; and one would expect 
this to lead to degeneracy of the solution. That is, these 
specific boundary conditions are expected to restrict the 
shape of the boundary itself, and thus in turn to restrict 
the topology of possible equilibria. Although the bound- 
ary conditions are specified, and the problem is thus rig- 
orously formulated, the above statement regarding the 
over-constrained nature of our problem nevertheless has 
not been shown to be useful in constructing a formal 
mathematical proof concerning the geometry of the con- 
figuration in the presence of a zero-line. Formally, it is 
easier to discard the TN for a given topology by direct 
construction of a solution with needed properties. Gen- 
erally, it is not clear at all how to construct a formal 
proof that the only solution of ([!]) with boundary con- 
dition Bj^ = for the type A topology is unique, and 
axisymmctric, thus defining the shape of the boundary 
itself! It is even less clear how to prove that there is no 
smooth solution for the type B topology, assuming that 
this statement is true. 

Finally, we note that one can simply produce artifi- 
cial discontinuities, but these are irrelevant to our dis- 
cussion. To illustrate, suppose that we place supercon- 
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ductive walls at the locations of the thick lines in both 
panels of Fig. 4. In that case, the configurations will be at 
equilibrium (for both type A and B topologies), and the 
field will be smooth everywhere except at the boundaries 
(where the field jumps from a finite value just outside the 
walls to zero at the walls in order to meet the boundary 
condition of zero field within the superconducting walls). 
This type of discontinuity is irrelevant to the astrophys- 
ical problem we are aiming at, and we therefore do not 
discuss it any further. 



B. Description of the solution method. 

One of the powerful ways to study the formation of 
current sheets is via numerical simulation. However, in 
numerical simulations the Lundquist number, S — caL/t) 
{ca the Alfvcn speed, and L the characteristic length), 
which is critical for this problem, is far below that corre- 
sponding to values encountered in natural systems, viz., 
under astrophysical conditions ||. When S is not suf- 
ficiently large, the separation between typical reconnec- 
tion times and typical fluid dynamical times may not be 
large; it is therefore difficult interpret realistic resistive 
calculations in the context of a problem in which current 
sheet formation is to occur without topological changes. 
On the other hand, numerical schemes which attempt 
to circumvent this problem by solving the ideal MHD 
equations suffer from the difficulty that such schemes 
may be subject to numerical instability, so that it be- 
comes difficult to distinguish between numerical artifact 
and physically correct current sheet formation. When 
discontinuities in the magnetic field appear, traditional 
numerical MHD codes tend to either break down, or to 
introduce a small amount of resistivity to broaden the 
current sheets (so that, for example, their width is larger 
than a mesh cell). In certain situations, the symmetry 
of the problem can be exploited to study the approach 
to the ideal solution, i.e., 77 — > 0, see, e.g., (m). How- 
ever, typical simulations actually add some amount of 
numerical resistivity, as in fl(i[| , so that numerical solu- 
tions of the ideal induction equation correspond to solu- 
tions of that equation with an added effective diffusivity. 
In studies of reconnection it is known that the specifi- 
cation of boundary conditions on the magnetic field and 
velocity (which specify the rate at which magnetic field 
and plasma is brought into the reconnection layer) may 
affect the rate of reconnection. In our simulations we 
study spontaneous formation of singularities by isolating 
the flux system from the boundaries. We surround our 
flux bundle by a (transverse) field-free region, and we 
place the boundaries far away from the bundle, thereby 
minimizing the effect of the boundary conditions on the 
formation of the current sheets. 

We address this issue as a relaxation problem in the 
framework of ICV flows. Our approach involves a direct 



numerical simulation of ICV flows, i.e., solving the set of 
equations ([}]) and the momentum equation, 

dw <9v , _. 

= --Vp+ — {V x B} x B + jA7 2 v, (4) 
p Anp 

with V • v = pi]. We use a Lagrangian approach to 
solve the induction equation as in 0, §. More 
specifically, the magnetic field inside the region of in- 
terest is represented by a large number of field lines; the 
evolution of the field lines is then followed using the exact 
Lundquist solution, i.e., knowing the initial strength of 
magnetic field on a fluid element connecting two nearby 
points on a field line, the final strength is proportional 
to the length of the segment, as it is stretched by the 
motions. We assume for all cases that the magnetic field 
vanishes on the outermost field line (the dashed curves 
shown in the figures). The number of field lines which 
fill the domain is chosen so that the subsequent field evo- 
lution can be followed without leaving gaps in the final 
state, i.e., we determine the number of initial field lines by 
fixing the spatial resolution of the final state; we discuss 
this point further immediately below. As an important 
aside, we note that the initial magnetic field is smooth, 
implying that the current system, y), which is defined 
by Ampere's law 

„ „ 47T . 

V x B = — j, 

c 

is smooth as well, i.e., there are no current sheets initially. 

The momentum equation (Q), in contrast, is solved us- 
ing standard finite difference techniques, with finite vis- 
cosity. However, the requirement of coupling the mag- 
netic field evolution to the momentum equation does lead 
to a complication for computing the Lorentz force. The 
key issue is that the momentum equation requires the 
Lorentz force to be evaluated on a homogeneous spatial 
grid, while the magnetic field evolution is given in La- 
grangian space. We resolve this issue by (quadratically) 
interpolating the Lorentz force at each time step onto the 
homogeneous grid used by the momentum equation (Q); 
similarly, we use quadratic interpolation from the mo- 
mentum equation mesh to evaluate the velocity field on 
the Lagrangian mesh. In order to minimize interpolation 
errors, we fix the number of field lines such that every 
Eulcrian grid domain is pierced by at least a few field 
lines throughout the calculation. Note here that inter- 
polation errors do lead to inaccuracies in the solution of 
the flow and magnetic fields, but by construction cannot 
lead to changes in the magnetic field topology. Note also 
that we have checked for convergence of the solutions as 
the spatial resolution of our calculation is increased; our 
conclusion is that the results presented here do not de- 
pend on grid resolution. Our solution corresponds to the 
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limit r\ — * 0, in the sense that the topology is strictly 
conserved, but with finite viscosity; thus, the computa- 
tional scheme we use forces the relaxation to be due solely 
to viscous damping, and as a consequence, the field re- 
laxes to an equilibrium state. Our approach has the dual 
virtues that the boundary conditions for the magnetic 
field do not need to be specified, and that the field topol- 
ogy is preserved; it is therefore appropriate for the study 
of TN. 

If the viscosity is large, then (||j4|) describe monotonic 
relaxation to equilibrium. We can estimate the relaxation 
time as follows: from (||) we find that v ~ c\L/v = caS u , 
where S v — caL/v. This viscous regime is realized if 
S y <C 1. The relaxation time is then t u ~ L/v = ta/S v , 
with t a = L/ca- In the opposite limiting case, S v 3> 1, 
the system undergoes (strong) Alfven oscillations (v 
ca), with a period ta, decaying on a viscous time t v ~ 
L 2 jv = taS v . These two cases can be jointly described 
by an interpolation formula, 

tu = T A {S U + l/Su), 

from which it follows that the relaxation time is large for 
both limiting cases (in terms of ta)- Thus, optimal relax- 
ation to equilibrium occurs for S v ~ 0(1); in the simula- 
tions, we used the value S v = 5. It is important that S v 
not be too large: An important constraint on the value of 
S u is that the simulations remain stable. This constraint 
is not met if S v is too large; because the two dynamical 
equations are solved in different coordinates [eq. (^|) in 
Lagrangian, and eq. (^|) in Eulerian coordinates], errors 
arise from the interpolation from one coordinate system 
to the other, and therefore the calculations make sense 
only if these errors are damped sufficiently by viscosity. 

III. DESCRIPTION OF THE RESULTS. 

We conducted two series of numerical experiments for 
case A. In the first series, we consider the relaxation of 
this type of topology without any external field, as de- 
picted in Fig. 1. We explored different initial shapes of 
the field lines, including ellipse-like, diamond-like, and 
other similar configurations. In addition, for a fixed 
shape, we explored different distributions of the flux func- 
tion A(x,y), i.e., different functional dependences A(s), 
where s labels the field lines. The results are always the 
same: the field ends up in an axially symmetric state, 
provided the field vanishes on the outermost field line. 

In another sets of experiments, this same configuration 
(case A) is placed between magnetic walls, as in Fig. 2a; 
in this case, the system always evolves to create disconti- 
nuities, as in Fig. 2b, where the field lines are taken from 
one of our simulation runs. We see that as the "bub- 
ble" evolves, it attempts to become axisymmetric, but 
as it does so, two discontinuities begin to form, as de- 
picted in Fig. 2b, see also Fig. 5a. It is interesting to 



note that, for some initial conditions, a current point is 
formed, rather than a current sheet (or a line in two di- 
mensions), suggesting that finite conductivity could pre- 
sumably result in fast reconnection; that is, according 
to the Sweet-Parker mechanism (see, e.g., ||), the re- 
connection rate va ~ where I is the length of the 
current sheet, so that a short current sheet speeds up the 
reconnection. (In the classical Sweet-Parker mechanism, 
t = L, and v d ^ca/S 1 ' 2 .) 




FIG. 5. Profiles of the z-current corresponding to (a) the 
topology of Fig. 2b, and to (b) Fig. 3b (the field line config- 
urations of Figs. 2b and 3b are reproduced at the bottoms 
of (a) and (b), respectively). The current sheets are repre- 
sented by negative currents. Due to the presence of an exter- 
nal zero line, the total current is zero, and therefore strong 
and peaked negative current is compensated by a spatially 
distributed positive current. Note that in panel (a), the cur- 
rent corresponding to the external field of the magnetic walls 
changes sign, because each wall contains two zero-lines. 

It is crucial to note here that the evolution we just de- 
scribed is not forced by the walls; thus, the field and fluid 
near the walls (i.e., on the wall side of the zero-lines) has 
the equilibrium property (||). To see this, note that on 
the zero-line, the total pressure is continuous. Thus, we 
could replace the initial elliptical field configuration (the 
"bubble" ) lying between the two zero-lines with a field- 
free region whose gas pressure exactly balances the total 
pressure on the wall side of the zero-lines. The resulting 
configuration is clearly in equilibrium, and makes clear 
that the walls do not push the bubble, i.e., that the evo- 
lution of the bubble is entirely driven by the fact that it 
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is not in equilibrium. Thus, it is as the bubble tries to 
become axisymmetric, and pushes back the walls, that 
the two discontinuities are formed. In principle, if the 
bubble could reach equilibrium with ellipse-shaped field 
lines as in Fig. 2(a), then it would not even interact with 
the walls, and the equilibrium of the whole configuration 
would be smooth. The field evolution to TN here de- 
scribed, i.e., evolution from an initially smooth state to 
a state containing a singularity, is therefore an intrinsic 
property of the initially smooth state, rather than being 
forced by external means. 

Consider now the type B configurations. We again 
conducted two series of experiments. In the first series 
of numerical experiments, we studied different kinds of 
initial states, with different initial field line shapes, and 
with different distributions A(s); in all cases, we again 
required that the outermost line must be a zero-line, as 
shown in Fig. 3. We found for the type B configuration 
that a field discontinuity always appeared, as in Fig. 3b 
(which is taken from one of our simulations), no mat- 
ter what the initial distribution of A, or what kind of 
analytical representation of the initial field lines we used. 

In the second set of experiments, we simulated the evo- 
lution of magnetic field with different number of field 
lines. The issue is as follows: The magnetic field gradi- 
ent at x = increases during the evolution, so that the 
current, V x B, approaches a ^-function (Fig. 5). It is not 
possible to observe this tangential discontinuity because 
in the simulations, the field is described via a finite (al- 
beit a very large) number of field lines (recall p8|). Our 
hypothesis is that in the limit of an infinite number of 
field lines N — > oo, the current at x = tends to infinity; 
in order to test this hypothesis, we increased the number 
N (recall Q) in a succession of simulations that were 
otherwise identical. According to our hypothesis, we ex- 
pected the current to grow roughly as 1/A, where A is 
the closest distance between the X-point and the nearest 
field line; the experiments confirmed this expectation. 



in which one expects to encounter small but finite resis- 
tivity, these flux systems are initially unlinked, but as 
they are pushed together (and begin to reconnect), flux 
linkage is expected to occur and to lead to a field topol- 
ogy analogous to that depicted in Fig. 2, or to the generic 
type B configuration, discussed here. The magnetic flux 
surrounding these two islands would be initially weak, 
and the current sheet which is formed is therefore ex- 
pected to be weak. However, during the course of recon- 
nection, more flux will be pushed outside the two islands, 
thus accelerating the process of reconnection. This pro- 
cess may therefore be self-accelerating, resulting in final 
(spontaneous) reconnection; preliminary numerical sim- 
ulations of a resistive case of this sort suggest that the 
reconnection rate Vd scales as ca/S 01 , where a is a small 
power, a ~ O(O.l) ^J. If confirmed, it would imply that 
the reconnection is fast enough to satisfy the observed 
(solar) constraints on reconnection times. (Recall that 
while the Sweet-Parker reconnection time for typical pa- 
rameters corresponding to the solar corona is about three 
years, the time corresponding to Vd — ca/S ' 1 is only 
30 minutes, which is comparable to the energy release 
time scale for large solar flares, related to the so-called 
"long-enduring" events S). Therefore the two topolo- 
gies depicted in Figs. 2 and 3 may be regarded as generic 
examples of "fast" reconnection and activity in magnet- 
ically active astrophysical systems. 
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IV. CONCLUSION. 



The fundamental result emerging from our simulations 
is that the vanishing magnetic field on the outermost field 
line imposes strict constraints on the geometry of equi- 
librium: The type A topology can be at equilibrium only 
if it is axisymmetric; and therefore, if constrained by ex- 
ternal walls, it is at TN. Similarly, the type B rosette 
structure develops discontinuities, but only in the pres- 
ence of an external zero-line. The presence of zero-lines 
is thus an important aspect of topological nonequilibria. 

Finally, we comment briefly on the applicability of 
these results to astrophysical situations. Observations of 
the solar atmosphere p(J commonly show topologically 
unconnected magnetic flux systems which are seen to in- 
teract (viz., emerging flux loops). In such circumstances, 
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